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Abstract 

Quantum wires and electromagnetic waveguides possess common features since their physics 
is described by the same wave equation. We exploit this analogy to investigate experimentally 
with microwave waveguides and theoretically with the help of an effective potential approach the 
occurrence of bound states in sharply bent quantum wires. In particular, we compute the bound 
states, study the features of the transition from a bound to an unbound state caused by the variation 
of the bending angle and determine the critical bending angles at which such a transition takes 
place. The predictions are confirmed by calculations based on a conventional numerical method 
as well as experimental measurements of the spectra and electric field intensity distributions of 
electromagnetic waveguides. 
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I. INTRODUCTION 



The opinion that a binding potential is necessary to obtain a bound state of a particle 



in quantum mechanics [l] was already exposed as being incomplete in the ear 
quantum mechanics. Indeed, in 1929 John von Neumann and Eugene P. Wigner 
that even a strongly repulsive potential can support a bound state with a normalizable 
wave function. Since then many examples of unusual bound states in various quantum and 
classical analogue systems with different bindung mechanisms have been investigated. For 
an overview see for example Ref. 

In quantum mechanics it is by now well-known that a binding potential in two space 
dimensions always supports a bound state no matter how shallow the potential, whereas 
in three space dimensions the potential needs to have a certain depth. The reason is that 
in two space dimensions the state of vanishing angular momentum, that is, a rotationally 
symmetric state, brings in an additional attractive potential which decays quadratically with 



y days of 
2| showed 



the inverse of the distance. This anti-centrifugal potential {4] arises from the adaption of 
the Laplacian to the symmetry of the state and appears in many other situations where the 
boundary conditions or the preparation mechanism dictate the symmetry of the state. In the 
present paper we develop an approximate but analytical technique to describe the occurrence 
of such geometry-induced bound states for the case of a bent waveguide. We compare and 
contrast our analysis with numerical calculations based on a conventional method as well as 
with measurements of these bound states. 



A. Prom compound nucleus to bent waveguides 

Systems of the type considered in the present paper in fact serve as models for the 
understanding of properties of the scattering processes associated with compound nucleus 
reactions. Within such a scattering approach the compound nucleus is modeled by a confined 
region of quasibound states, and a number of open reaction channels attached to it. These 
represent the coupling of the quasibound states to the continuum, i.e., their formation and 



decay jsNlol [see Fig. l(a)| . 



The universal scattering properties of a compound nucleus reaction, in turn, can be sim- 
ulated in microwave experiments, where electromagnetic waves propagate via antennas or 
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waveguides into a cavity corresponding to the interior region which allows for bound or qua- 



sibound states 



lOHlSj. The most elementary example is a bent waveguide with a corner as 



depicted in Fig^ 1(b) — a model for strong interactions in nuclei similar to that one was pro- 
posed in Ref. jl4|. States trapped around the inner corner were predicted to exist below the 



energy of the first propagating mode 14j-|l8[| and o 



waveguides jl9l-l2l| and in quantum wires 
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3served experimentally in bent microwave 



23 



26| . Note that the Helmholtz equation 



describing the physics of the former and the Schrodinger equation which is applicable to the 
latter are mathematically identical in the energy range of the bound states. This analogy 
was used in the numerical computation of bound states in bent waveguides j20] and quantum 



wires 
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26[. The experimental and the numerical studies revealed that with the variation 
of the waveguide's bending angle a [see Fig. l(b)| a transition of the modes from bound 
to propagating ones takes place. The bending angles at which such a transition occurs are 
called critical angles in the following. 



B. Formulation of the problem 



In this paper we present a detailed theoretical and experimental study of that transition. 
To test the analytical approach, existing results [20| were extended down to bending angles 
a = 0.2° using an improved conventional numerical method which is based on a scattering 
approach. Its main achievement is an analytical formula for the number of bound states at a 
given bending angle derived from a fit to the numerical result. The experiments performed to 
thoroughly test these results comprised measurements of the eigenfrequencies (energies) and 
the associated field intensity distributions (squared wave functions) of the bound states for 
several bending angles. Furthermore, we present experimental results concerning the effect of 
the finite lead lengths on the bound states and the occurrence of resonant tunneling. These 
investigations go beyond the numerical and theoretical ones which assume waveguides of 



infinite lengths. Only recently the existence o 
in the experiments, was studied theoretically 



trapped modes in finite waveguides, as used 



221 



The numerical computations of the bound states and the critical angles of a waveguide 
or quantum wire involve the solution of the Schrodinger equation of a free particle with 



Dirichlet 
17 



joundar y co nditions imposed on its wave functions at the walls of the system 



18, 



20 



15 



24, 



25, 



28 



29|. 
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The theoretical treatment of a problem of this kind presented in Ref. 30|] encompasses 
a conformal mapping applied to the geometry of the bent waveguide which induces in the 
Schrodinger equation an additional effective potential given by the associated Jacobian. 
Starting from this equation the bound states at a given bending angle and the critical 
angles were determined using a Wentzel-Kramers-Brillouin (WKB) approximation developed 



in Ref. 



31| for systems containing sharp corners. The results agree well with those obtained 



experimentally with electromagnetic waveguides and from the numerical computations based 
on a scattering approach. 

C. Relation to other work 



It is worthwhile to mention again that some of the salient features of the spectrum of 
a bent waveguide were already understood in previous experiments with electromagnetic 



waveguides [19142 1| and quantum wires [22|, l23|. Nevertheless, we feel that it is important 



to obtain a complete solution, by improving the conventional numerical and experimental 
methods, which encompasses all the important features and providing a general method 
which does not rely on numerical calculations. For example, one may pose the following 
question: What is the number of bound states supported by an orthogonal cross-wire con- 



figuration? The numerical answer is two 



15j . but the reason for this should be revealed by 



a suitable treatment relating the shape of a wire to its binding capabilities. 

These effects, rather than being just accidental, are linked to deeper concepts in undu- 
latory physics. The confinement of a particle to a restricted domain by introducing walls 
is directl y re lated to a proper formulation of the d'Alembert principle in quantum mecha- 



nics 
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33|. This problem must be carefully tackled in the light o: 



such as the trapping of waves in smoothly bent quantum wires 



emergence of a so-called curvature-induced quantum potential 



It oi weil- 



1-known results 



34| and the 



35| . For instance, it is strik- 



ing that in bent waveguides or quantum wires with open ends as depicted schematically in 
semiclassical methods such as trace formulas derived for closed 0, Q and open 



Fig. 1(b) 



systems 
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B|-|40| fail because there are no closed orbits and there is no classical counterpart 
including a binding force. We present an analytical approach which yields a Hamiltonian 
whose spectrum comprises the bound states. The conclusion to be drawn from this anal- 
ysis is that the mere presence of curved boundaries may induce a binding potential, with 
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significant implications close to sharply bent boundaries. 



D. Outline of the article 

The article is organized as follows: in Sec. |TT] we shortly review the results of our cal- 
culations of the bound states in a bent waveguide and of the critical angles based on a 
conventional method. Section IIIII comprises the description of the microwave experiments 
and the comparison of the experimental and the numerical results from Sec. [TTl Section IIVI 
presents analytical estimates of the critical angles and the number of bound states in a bent 
waveguide. The experimental and the numerical results presented in Sees. HTl and IIIII demon- 
strate that for a given bending angle the wave functions of the bound states with excitation 
frequency close to the cut-off frequency penetrate considerably into the leads. Similarly, the 
influence of the binding potential revealed in the theoretical approach of Sec. IIVI reaches 
into the leads. In Sec. |V] this phenomenon is investigated experimentally. The conclusions 
are summarized in Sec. I VI I 



II. NUMERICAL RESULTS 



In this section we compute the eigenenergies and wave functions for the bound states of 
sharply bent quantum wires, that is, the eigenfrequencies and the associated electric field 
strengths for those of sharply bent rectangular waveguides using a conventional method. 
For this we use the fact that below the frequency /max = c/{2h), with c the speed of 
light and h the height of the waveguide, only the TEi o-mode is excited 4l|. There, the 
Helmholtz equation associated with a bent waveguide is mathematically identical to the 
stationary Schrodinger equation of an open quantum billiard or a quantum wire of corre- 



sponding shape 



11 



42 



43| . The eigenfrequencies of the bound states are all below the cut-off 



frequency fc = c/{2W) of the first propagating mode, where W is the width of the wave- 
guide. Consequently, since rectangular waveguides with the width larger than the height 



were used in the experiments, such that /max > fc, the numerical solutions of t 



Schrodinger equation obtained in this section apply to bent quantum wires 
to bent microwave waveguides. 



le stationary 



23 



26| and also 



The stationary Schrodinger equation of an open quantum billiard with the shape of a 
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bent waveguide corresponds to that for a free particle with Dirichlet boundary conditions 



along the boundary dil defined by the walls of the waveguide in Fig. 1(b) Here, denotes 
the interior of the waveguide. A coordinate system (x, y) is introduced in the plane of the 
waveguide with origin at the outer corner and the x-axis along the symmetry axis of the 
waveguide to write for x, ?/ G 

- A^^y(j){x,y) = E(f){x,y), = 0, (1) 

where we use units fi^ / (2m) = 1, the energy E corresponds to the square of the wave number 
k, E = k'^, the wave function 0(a;, y) denotes the 2;-component of the electric field vector in 
the waveguide and 



is the two-dimensional Laplacian. The width of the waveguide is set to ly = 1 in all 
calculations. 

The procedure for the numerical determination of the eigenenergies and eigenfunctions 
of Eg. ([Tl) is outlined in Appendix |Xl It is based on the scattering approach described in 



Refs. 



24 



25 
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44l |. For this purpose the waveguide is subdivided into an interior region 



and an asymptotic region defined by its straight leads with parallel walls as indicated in 



Fig. 1(b) The procedure is similar to that used for the description of compound nucleus 
reactions in the framework of scattering theory, where the interior region corresponds to the 
compound nucleus and the leads to the formation and decay channels jo], . 

Figure [2] shows the numerical values k/k^ = f/fc for the bound states as function of the 
bending angle a. Here, k = \/E is the wave number, and / = ck/ {27t) with c the speed 
of light is the corresponding frequency. At the cut-off frequency fc = ckc/ {2tt) = c/ {2W) 
the first propagating mode emerges in the leads [see Eq. ( ]A4[) ]. Each curve corresponds to 
the eigenfrequency of one bound state which varies continuously with the bending angle a. 
The number of bound states at a bending angle a is obtained by determining in Fig. |2] the 
number of curves at a. All curves start at a = 0° and increase until they eventually reach 
the value ///c = 1 at a critical angle acrit- There, the number of bound states decreases by 
one. Above a^^it ~ 28.28°, with denoting the critical angle of the nth bound state, only 
one bound state survives. It corresponds to the lowest curve in Fig. [2] and has a^j^ = 180°. 
With decreasing a the number of bound states at a given value of a increases rapidly. 



Indeed, as can be observed in Fig. 2(a), it is hardly possible to resolve all curves around 



a ~ 0.2° due to their large number. For a — )■ 0° the eigenfrequency of the lowest bound 
state approaches from above the value ///c = 1/2, so all eigenfrequencies obey the inequality 
f I fc > 1/2- Note that the frequency / = /c/2 corresponds to the cut-off frequency of a 
straight waveguide with twice the width of the leads. Indeed, such a waveguide is obtained 



by decreasing the angle of a to 0° in Fig. 1(b) 



Figure [3] shows the intensity distributions corresponding to the four bound states existing 



at a = 9.65°. These are marked by dots in Fig. 2(b) For the first bound state at ///c = 0.631 
the intensity distribution shown in panel (a) exhibits one maximum close to the inner corner. 
With increasing ///c, that is, for the second [panel (b)], the third [panel (c)] and the fourth 
[panel (d)] bound state the number of maxima increases by, respectively, one. The maximum 
with the highest intensity moves away from the inner corner to the outer one, and that 
closest to the inner corner penetrates more and more into the leads. This leakage has been 
investigated experimentally. The results are presented in Sec. |Vl 

With the numerical method presented in Appendix |A] we were able to compute the eigen- 
frequencies of the bound states of bent waveguides with bending angles as small as 0.2°. For 
smaller angles the numerical errors are larger than the spacing between the eigenfrequencies 
with values / ~ /c. 

In Fig. mthe number of bound states A/'(a) is plotted versus the bending angle a. The 
circles correspond to the numerical results. They are well described by the analytical equa- 
tion 



(3) 



Ar(a) = lnt[ycsc (f)+V2 

with the fit parameter oq = 1.027 (solid line) close to unity and Int[x] denoting the integer 
part of X. 



III. MEASUREMENT OF BOUND STATES 



In the preceding section we have presented numerically determined eigenfrequencies and 
associated wave functions, i.e., electric field strengths, of the bound states in bent rectangular 
waveguides as function of the bending angle and determined the critical angles, where a 
bound states turns into a propagating mode. These results were tested experimentally with 
bent microwave waveguides as described in the present section. 
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A. Experimental setup 



Experiments with microwave cavities are widely used to investigate universal properties 



of quantura systems such as the spectral statistics of chaotic quantum billiards |46|, chaotic 
scattering |lol, or fidelity decay [48;]. This approach is possible due to the analogy noted 
already above between the stationary Schrodinger equation and the Helmholtz equation for 



fiat cylindrical microwave cavities 



11 



42| . Similarly, the properties of quantum wires have 



3een investigated in experiments with microwave waveguides of rectangular cross section 



19 



-|21|. Here we report on microwave experiments performed to validate the numerical and 
the theoretical results concerning the critical angles and the dependence of the bound states 
on the bending angle. 

Nine different waveguides with different bending angles a were investigated experimen- 



tally. Figure 1(b) shows the geometry of the bent waveguides. They were constructed by 
soldering together two tapered WG18 brass waveguides (by Flann Microwave) as shown in 
Fig. [51 The leads with lengths /q were extended to their full lengths / by attaching additional 
straight waveguides with open ends. The angles of three of the nine waveguides were chosen 
close to a critical angle, a^^j^ ~ 28.28°, ol~^)^^ ~ 16.26° and a^^it ~ 11.48°, which were read 



off from Fig. 2(b) and are marked in this figure by arrows. The lengths / were chosen larger 
than 400 mm, in order to minimize the effect of the finite length of the waveguides analyzed 
in more detail in Sec. |Vl The angles and the lengths of the different waveguides are listed in 
Tab. [B They were determined with a precision Aa = 0.1° and AZ = 0.5 mm, respectively. 
The waveguides had a rectangular cross section with a width W = (15.76 ± 0.02) mm and 
a height h = (7.9 ± 0.02) mm. 

According to Eq. flA4l) . a straight waveguide of this width allows for propagating modes 
above the cut-off frequency fc = (9.511 ± 0.012) GHz. Below fc, only evanescent fields can 
exist. Note that the quantities of interest, the ratios f / fc, do not depend on the width W 
of the waveguide. Hence, the experimental ratios may be compared with the computed ones 
of Sec. HI] for all choices of W. 

Below /max = c/ (2/i) ^ 19 GHz the electric field is perpendicular to the waveguide plane 



defined by the x and y coordinates and described by the scalar Helmholtz equation j43|, |49 | 

{A,,y + e)E,{x,y) = 0, E.Ian = 0. (4) 
Equation (HD is identical to Eq. (II]) with (p = and E = = {27if/cf. 



Microwave power was coupled into and out of the waveguides with two small wire antennas 
aligned perpendicularly to the waveguide plane. They protruded about 1.5 mm into the 
waveguide through holes with a diameter of 3 mm in its top wall at three different positions, 
one in the interior region (label 1 in Fig. E]) and two close to it in the leads (labels 2 and 
3 in Fig. E]). Results are presented for the measurements with antennas at positions 2 and 
3. A vectorial network analyzer (VNA, model PNA N5230A by Agilent Technologies) was 
used to measure the transmission amplitude Sba{f) from antenna a to antenna b, where 

\SUf)\' = ^ (5) 

is the ratio of the power Pout, b coupled out via antenna b and the power Pi^ ^ coupled in via 
antenna a at the excitation frequency /. 



B. Resonance spectra and frequencies and field intensities 

The measured resonance spectrum of the waveguide with bending angle a = 16.40° is 
shown in Fig. |6l There are three sharp resonances labeled by (a)-(c) below the cut-off 
frequency (gray bar) at frequencies 6.49, 8.16 and 9.50 GHz. Their quality factors Q ~ 2000 
are comparable to those of resonances in closed normal conducting microwave cavities of 
similar size, which indicates that these resonances correspond to the three bound states 
predicted for this waveguide. Above fc the spectrum exhibits a series of broad resonances 
like the one labeled by (d). These are due to partial reflections of the microwaves at the 
open ends of the waveguide. 

To obtain a better understanding of the resonances their electric field intensity distribu- 



tions were measured using the perturbation body method [50|,|5l|]- A cylindrical perturbation 
body with diameter 2 mm and height 2.6 mm made of a rubberlike plastic combined with 
magnetic barium ferrite powder was introduced into the cavity and moved by an external 
guiding magnet to different positions {x,y) on a grid with 2.5 mm spatial resolution. For 
each position of the perturbation body, the shift A/ of the resonance frequencies was mea- 



52|, / cx \Ez{x,y)\'^, at the position 



sured. It is proportional to the electric field intensity 
of the perturbation body, 

A/oc |E,(x,2/)|2. (6) 
The measured intensity distributions at the frequencies of the resonances with labels (a)- 
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(d) in Fig. [6] are plotted in Fig. [71 In panels (a)-(c) they exhibit yellow (bright) stripes, the 
color of which does not change to blue (dark) towards the walls of the waveguide, that is, 
the intensity does not seem to vanish there. This is an experimental artifact resulting from 
the friction between the perturbation body and the surface of the brass waveguides, which 
can lead to a lag of it behind the guiding magnet yielding an inaccuracy in its position. 

The intensity distributions shown in panels (a), (b) and (c) of Fig. [7] are localized in the 
corner region of the waveguide and decay in the leads of the waveguide. The resonances 
(a) and (b) evidently correspond to bound states as expected, since their frequencies are 
well below the cut-off frequency. The frequency of resonance (c) lies within the gray bar 
indicating the range of the expected value of fc (see inset of Fig. E]). Its intensity distribution, 
however, is localized in the corner region, which demonstrates that it also is a bound state. 
Still, it extends significantly into the leads. This is further discussed in Sec. |Vl Note that 
the length / of the leads was chosen such that it is much larger than the decay lengths of 
the electric field intensities associated with the bound states. 

In contrast to the intensity distributions associated with the resonances (a)-(c), that 
corresponding to the resonance (d) extends over the whole waveguide and shows a periodic 
behavior in the leads [see Fig. [3(d)]. This is expected since its resonance frequency is larger 
than the cut-off frequency fc for propagating modes inside the leads. The resonance itself is 
very broad with Q ^ 200 since microwave power leaks out at the open ends of the waveguide, 
i.e., it corresponds to an unbound state. Figures [Tl^a)-(d) demonstrate that bound and 
unbound states can be clearly distinguished by means of their intensity distributions. 

Figure [S shows the measured eigenfrequencies in units of the cut-off frequency fc of the 
bound states (crosses) together with the values calculated in Sec. [Ill with a conventional 
method (solid lines, cf. Fig. [2|). The agreement of both is good except for small deviations 
which are attributed to imperfections in the construction of the waveguides resulting from 
the soldering process. 

It should be noted that all the measured frequencies are below fc, even though the 
bending angles a = 28.40°, 16.40° and 11.83° are slightly above the corresponding critical 



angles marked by the arrows in Fig. 2(b) However, in all three cases the frequency of 
the highest bound state is within the error range Afc = 12 MHz of the cut-off frequency 
indicated by the gray bar in Fig. [6l and the bending angle of the waveguides is only known 
with a precision of Aa = 0.1°. In conclusion, within these uncertainties the measurements 
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confirm the values of tlie critical angles deduced from the numerical calculations of Sec. HIl 



IV. EFFECTIVE POTENTIAL FOR SHARPLY BENT WAVEGUIDES 

In this section we present an analytical approach to determine the bound states in the 
waveguides and the values of acrit- It involves the solution of the Schrodinger equation for 
a Hamiltonian with an effective binding potential, induced by a conformal map, based on 



an extension of the WKB method to systems comprising sharp corners 3l|. A first work 



on this subject was published in Ref. [30|, which unfortunately contains misprints. For 
this reason we revisit the approach and derive an improved approximation for the effective 
binding potential. 

A. The conformal map 

The essential idea of the theoretical approach in Ref. [s^ was to transform the Schrodinger 
equation ([T]) with the coordinates x,y & Q confined to the interior of a sharply bent wave- 
guide to one with coordinates < m < 1, — oo < f < oo restricted to an infinite strip. 
One procedure to obtain such a transformation is to introduce complex variables z = x + iy 
and ^ = u + iv and a transformation z = F{C,), which maps the infinite straight strip in 
the complex plane defined by the range of values of u and v onto the interior of the bent 
waveguide in the (x, y) coordinate system. 

There are many possibilities to define F{^). However, we impose the requirements that 
it should be smooth everywhere except at the corners of the bent waveguide and that the 
Laplacian A^^y = + ^ should be transformed into a Laplacian A„_„ = ^ + ^ without 
cross terms. 



For a > 0° a transformation with these properties was constructed in Ref. 30|| using a 
combination of two Schwarz-Christoffel mappings [53i], 



e^X = sin2(^e) ^^ = i^(0 (7) 



which images the infinite strip ^ = u + iv onto 

^ = m = ^/o'^dttt-i(i-t) 
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where a = and B{x,p,q) is the incomplete Beta function |54|. 

The image z{u,v) = F{C,) of the point {u,v) = (0,0) is the outer corner of the bent 
waveguide at {x,y) = (0,0), that of the point {u,v) = (1,0) is its inner corner at {x,y) = 
(csc[a/2], 0), while the line (0,f) is mapped onto its outer wall and the line (l,f) onto its 
inner wall. 

Figure [9] shows three bent waveguides with angles a = 180°, 126° and 54°, together 
with the coordinate grid {x,y) = {x{u,v),y{u,v)). For a fixed v = const, the coordinate 
lines approach the boundaries along {x{0,v),y{0,v)) and {x{l,v),y{l,v)) perpendicularly, 
whereas for a fixed u = const, they are parallel to the boundaries for large values of v. In 
Fig. [9] examples are shown for both types of coordinate lines. 



B. The Jacobian 



The change of coordinates implies the transformation 

d{u,v) . 



A 



where the inverse of the Jacobian 



""'^ d{x,y) 
d{x,y) 



d{u, v) 



dz 



appears as a prefactor and the Schrodinger equation ([T]) takes the form 

d{u, v) 



d{x,y) 

with ip(u, v) = (f>[x{u, v),y{u, v)] and 



Au,vip{u,v) = Eip{u,v) 



1^(0, v) = ip{l,v) = 0. 



(9) 



(10) 



(11) 



(12) 



To compute the Jacobian we perform the variable transformation t = sin^ (f^^) 
integral of Eq. (|8]) and obtain the expression 



Jo 



cot ( -^9 



l-a 



which with Eq. ( iTOl) yields 



J{u, v; a) 



cot 



T 1-Q 



' cosh(7r?j)+cos(7r?i) 
^ cosh(7rj;)— cos(7rji) 



1-Q 



(13) 

(14) 
(15) 
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For (5 = 1, i.e., a = 180° we find J'{u,v;a = 1) = 1 and the Schrodinger equation reduces 
to that of a straight waveguide, the solutions of which are known. 

The abrupt change of the contour density observed in Fig. M for a < 180° close to the 
corners is related to the behavior of the Jacobian in their vicinity. In fact, near the outer 
corner, where |^| ~ 0, the Jacobian can be approximated as 

J{u,v;a) ~ (^j ' (16) 

that is, it has a singularity of the type a;""^ with the order 0<7 = 2 — 2q:<2. On 
approaching the inner corner, i.e., for |^ — 1| — t- 0, 

2 



J{u,v;a)^(%-^\ (17) 



becomes vanishingly small. For large values of f ^ 1 the Jacobian can be approximated as 

J{u, f ; 5) ~ 1 + 2 (1 - a) e'"''' cos(7rM) (18) 

and converges to unity. This behavior is expected since there the transformation from ^ to 
z = F{^) simply corresponds to a rotation of the coordinate grids by an angle | (1 — a) 
with respect to each other around the inner corner. 



C. Emergence of the potential 

Multiplication of the Schrodinger equation ffTTl) with the Jacobian J'{u,v; a) yields 

(-A,,„ + [1 - J{u, v; a)] E) ^{u, v) = Ei;{u, v), (19) 

that is, the mapping Eq. ([7]) entails the emergence of an effective binding potential depending 
on the Jacobian J'{u,v;a), i.e., the area density. To obtain an approximate analytical 
solution of this boundary value problem it is turned into an effective potential problem 
in the v coordinate by choosing an appropriate orthonormal basis for the u coordinate. 
The M-coordinate lines approach the f-coordinate lines at the boundaries of the waveguide 
perpendicularly (see Fig. [9]). Thus a suitable choice that obeys the boundary conditions 
imposed on i/j{u,v) [see Eq. ( ITT]) ] is given by 

oo 

(v) sin (nvrw) (20) 

n=l 

13 



with 

fn{v) = y/2 / dw f) sin (nvrn) . (21) 
Jo 

Furthermore, since the only singularity of the Jacobian, which is located at {u, v) = 
(0, 0), has an order less than 2, products of the form sin(m7ru) sin{n7[u)J'{u, v; a) are regular 
functions at all points inside the waveguide. Therefore, we may expand the function ■?/'(«, v) 
in Eq. f|T9|) in terms of the orthonormal functions a/2 sin(?7,7ru), n = 1,2, . . . and obtain the 
system of coupled equations 

oo 

E [-^nmd', + EVnUv)] fm{v) = EJ^{v) (22) 

with 

En = E - {mrf , Vnmiv) = 5nm " Vnmiv), (23) 

where 

Vnmiv) = 2 / (\usm{nTiu)sm{mTiu)J{u,v]a) (24) 
Jo 

are the coupling matrix elements. In this equation the energy E appears parametri- 
cally as the strength of the coupling matrix elements Vnm{v) which define an infinite- 
dimensional matrix. The numerical evaluation of the integral Eq. ( 12^ yields three char- 
acteristic properties of Vnm{v) which may also be deduced from those of the conformal 
map: (i) the matrix elements Vnm{v) decrease as one moves away from the diagonal, 
|V^,„+m('i^)| ~ \Vnn{v)\e~'^'" , m = 1, 2, ■ ■ ■ for all v, {ii) the diagonal elements saturate with 
increasing n, i.e., Vnniv) < Jq <\uJ{u, v; a), and (iii) for t> — oo the matrix Vnm approaches 
the identity matrix, Vnmiv) — > 5nm, i-e., Vnmiv) 0. 

D. Effective potential 

Numerical tests revealed that for small values of v ^ 0.25 we may write to a good 
approximation 

n-l 

E iVnmiv) + Vn,niv)) + Kn(w) ~ uV^niv), (25) 
m=l 

where the exact range of validity depends on n and on a, whereas for f ^ 1.5 

Vnmiv) ^ Vnniv)5nm- (26) 
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Accordingly, for small and large values of v we may replace Eq. (122|) by 

[-dl + EVn{v)) Uv) = E„Mv). (27) 

with 

Vniv) = 1 - nVn{v) for V < 0.25, (28) 

Vn{v) = l-Vn{v)f0TV>1.5, (29) 

where we use the notation Vn{v) = \4„(f). 

In terms of a Hamiltonian H{v) the Schrodinger equation would read 

[H{v) - E^] Uv) = 0, (30) 



with the eigenenergies En defined in Eq. (1231) . They are related to the eigenvalues E of 
Eq. (ITT]) which appear in Eq. (1271) as the pref actor of the potential Vniv). 

The coupling matrix elements defined by Eq. fl2^ depend on the bending angle a. Fig- 
ures 



10(a) and 10(b) show Vn=i{v) = 1 — Vn=i{v) and its derivative —^Ya^J^ for six different 



values of a. The potential Vn{v) is a symmetric function of v, has a minimum at f = 
and approaches for l^;! ^ 1.5. Furthermore it has the shape of a kink for a ^ 30°. Sim- 
ilarly for decreasing a its derivative, though vanishing at f = 0, undergoes an increasingly 
rapid change from its minimum V^^^ = — to its maximum V^^^ = — ^^^^^ situated 

at t> = ±e with e ^ 1, which is comparable to a jump k = V^^^ — V^^^ for small bending 
angles a. 

E. Semiclassical approximation 
In Ref. 



3l| an extension of the WKB approximation to systems described by a potential 
which displays a corner and has a discontinuous derivative is proposed. We now apply this 
approach to the present situation and determine approximately the bound states. 

1. Quantization condition 

Due to the minimum of the effective potential, bound states may exist for £"„ < EVn{v) < 
0. The value of v, where equality En = EVn{v) holds, defines the turning point Vt = \v\. 
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For n = 1 the inequality is not satisfiable for E > tt'^. This value defines the cut-off energy 
Ec above which no bound states exist. For E = 0.99997r^ we obtain vt — 1.5 for n = 1. 
For n > 2 the value of the turning point is smaller than 0.25. Furthermore, the extended 
WKB approximation involves the potential around the corner at f = 0. Accordingly, we 
may replace Vn{v) by 1 — nVn{v) in Eq. using Eq. 

The corner-corrected WKB approximation yields for a symmetric potential the quanti- 



zation condition 



3l| 



A 

— = tan 
4 



(31) 



where we defined the action 



S = £ dv^{n'Ky - EnVn{v), (32) 



and introduced the notations 



A 



\2 



-3/2 K 



EnVn{^ -nE, (33) 

and A = 0,1,2,.... Since the solutions of Eq. (!3T]) do not depend on the value of A we may 
set A = 0. 

Although the explicit forms of the potential and its derivative, which is needed for the 
computation of k, are not available, the associated integrals Eq. f l24l) and thus Eq. ( l3Ti) can 
be solved numerically. For simplicity we define the rescaled energies e = ^^^j^^ and Eq. (155]) 
takes the form 

_ .^3/2 ^ _^ e cos(g) + sin(g) 

^ ' 2K(0)47ry^cos(5) -sin(5)' 

For n = 1, 2, . . . , A/fe(a), with A/'b(a) denoting the total number of bound states at a, 
this equation yields the rescaled eigenvalues e„(a) of the bound states of a bent waveguide. 
Here, the values of e„(a) increase with their index, i.e., e„(Q;) < e„+i(a). They are bounded 
from above by the condition e < A/f,(a). Beyond that value all states are unbounded. Thus, 
the cut-off value of the rescaled energies is given by eda) = A4(a)- 



2. Number of bound states 



We determined the number of bound states A/'b(a) with the help of [55|, l56| 



/ dudpu J dvdpv 



27r 



- du r dv^^. 

TT Jo J-oo 2 



(35) 
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The value of Pmax is obtained from the classical analog of Eq. ( !T9|) . which can be written as 
-Eciass = Y + (1 — J{u, v; a)) and satisfies the inequality -Eciass < for the bound states. 
Thus, = {J{u,v;a) - 1) ^, yielding 

A4(a) = - f^du r dv {J{u, v]a)-l) + l = + 1. (36) 

TT JO J-oo TT 

Note, that for finite lead lengths, i.e., for v restricted to some value v < t'maxj 
du Jl'^^^^^ dv {J{u, V] a) — 1) = A{a) — 2finax, where integration over the Jacobian, A{a), 
yields the area of the finite bent- waveguide. The integral already approaches the value A{a) 
for t>inax — 1-5. In fact, A{a) can be interpreted as the effective area covered by the potential, 
as it corresponds to the area of an unbounded bent waveguide minus the area of the leads 
where the potential becomes vanishingly small. In the interval considered, 0.7° < a < 30°, 
it is well described by A{a) = csc(a/2) — qq, with oq — 0.91, that is, to a good approximaton 
the integer part of A/b(a) coincides with J\f{a) given in Eq. ([3]). It equals for a = 180°. 



3. Approximate critical angles and rescaled energies 

For a < 180° the rescaled energies, i.e., the ratio of the solutions e„(a) of Eq. and 
the cut-off energy eda) is given by 

en(a) = -r7^e„(a). (37) 
A/ [a) 



Generally the dependence of the solution of Eq. flMl) on 5* is negligibly small for a > 1°, 
so the last term in this equation can be set to unity. Then the solutions of Eq. fl34p are given 
as the zeros of a polynomial of third order. These are computed with the help of Cardano's 
method js^]- 

To determine the critical angles, where a bound state turns into an unbound one, we 
furthermore have to take into account that Eq. flM|) does not yield the correct ground 
state energy [3l|. Indeed, for a — )■ 0° the rescaled eigenenergy 61(0;) = ei{a)/ec{(y) does 
not converge to the expected value e/cc = P / fl = 0.25 (see Fig. [2]). However, when 
introducing a constant shift Cq ~ 0.25 of the rescaled energies e„(a), which ensures the 
correct convergence for a — )■ 0°, a good agreement between the solutions Ocrit of the equation 
Co + ^n{oicrit) = 1 a^d the numerical results of Sec. [TTl and thus also the experimental ones 
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of Sec. mil is achieved. This property is also true for the bound states as demonstrated in 



Fig. 11(a) Aheady for the second bound state the agreement between the numerical values 
of f / fc of Sec. ini and those obtained with the corner-corrected and ground-state shifted 
WKB approximation is good. 

The relative deviations of the critical angles obtained with the two methods are shown 



in Fig. 11(b) They are largest for large critical angles, however less than 2%, and decrease 
rapidly to a relative deviation of less than 1% for ctcrit ~ 1°- Thus the WKB approximation 
proposed in Ref. jsil for systems with a potential displaying a kink provides a good prediction 
of the critical angles. However, it is necessary to include the known result for the behavior 
of the rescaled energies, respectively, frequencies in the limiting case a — )■ 0°. 



F. Summary 

In summary, the conventional numerical methods, also the one presented in Sec. HIl used 
to compute the eigenfrequencies and wave functions of the bound states involve the solution 
of the Schrodinger equation of a free particle with Dirichlet boundary conditions, Eq. ([1]). 
In contrast, in this section a Schrodinger equation for a Hamiltonian including an effective 
potential was derived, the spectrum of which contains the eigenfrequencies of the bound 
states in a bent waveguide. These are determined by solving the associated Schrodinger 
equation fl27p with a semiclassical approximation. A comparison with the numerical results 
of Sec. [TTl and the experimental ones of Sec. IIIII corroborates the validity of the theoretical 
approach in terms of an effective binding potential presented in this section. The property 
of the effective potential that it is nonvanishing for v ^ fmax implies that its influence 
reaches into the leads of the waveguide, in accordance with the penetration of the wave 
functions into the leads as observed in the calculated intensity distributions in Fig. |3] and 
the experimental ones in Fig. [71 This feature already indicates that the existence of a 
bound state should depend on the length of the leads of the waveguide. The results of an 
experimental investigation of the effect of the finite length are presented in the next section. 
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V. FINITE LEADS AND RESONANT TUNNELING: EXPERIMENTAL INVES- 
TIGATION 

The numerically calculated intensity distributions and the experimental ones in Figs. [3] 
and [3, respectively, clearly demonstrate that there is a leakage of the bound state wave 
functions from the interior region into the leads of the waveguide. This section focuses on 
the experimental investigation of the decay lengths of the bound state wave functions into 
the leads. 



A. Decay into the leads 



While most theoretical works on bound states in bent waveguides assume leads of infinite 
lengths, this oversimplification of course is not feasible in the experimental realizations. Due 
to their leakage into the leads, modes localized in the corner region of a bent quantum wire 



or waveguide of finite length may coup 
reservoir in the case of quantum wires 



e to a continuum of states, e.g., to a 2d electron 
or to the electromagnetic waves in 



22 



23 



26 



free space in the experiments presented in this paper. 

Similarly, waves sent into such a device may be trapped in the interior region before they 
are scattered back to the exterior. This phenomenon manifests itself in transmission spectra 



measured between both 
via the bound states 15 



ead ends as resonances and is attributed to resonant tunneling 



22 



23 



26 



58 



59| which can significantly alter the transmittance 



characteristics of the devices. Strictly speaking, in a waveguide coupled to a continuum 
the associated states are no longer bound but acquire a finite lifetime since they can decay 
via tunneling through it. In this section we show some qualitative results on the effect of 
a finite lead length and investigate the transmission through the waveguides by attaching 
waveguide-to-coaxial adapters to the ends of the waveguides. 

The measurements of the frequency spectra and the intensity distributions presented in 
Sec. mil were repeated with waveguides without extensions, i.e., with lead lengths Iq given in 
Tab. [land open ends (as in Sec. IIIip . The decay lengths of the field intensity distributions 
were determined experimentally from the projection of the measured intensity distributions 
onto the coordinate xl parallel to the waveguides [see Fig. l(b)| . A so-obtained curve is 
depicted in Fig. [T21 where the measured intensity distribution of the second bound state of 
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the waveguide with bending angle 27.65° at 9.364 GHz is shown. 

The projection of the intensity distribution onto xl is represented by a gray hne in hnear 
(top panel) and logarithmic scale (bottom panel). Three different regimes are visible in 
the logarithmic plot (bottom panel). The line is slightly curved for ^ 20 mm, it is 
linear for 20mm ^ ^ 100 mm, i.e., I-Ezp decays exponentially, I-Ezp oc exp {—Xl/6), in 
the leads of the waveguide and for xl ^ 100 mm the noise level is reached. The decay 
length was determined with a linear fit (dashed black line) in the middle part of the line as 
^expt = (16.2 ± 0.4) mm which agrees well with the value 6 = 14.4 mm computed from 

5 = j^—. (38) 

In this case, the eigenfrequency of the bound state is well below the cut-off frequency and 
the associated intensity distribution extends not very far into the waveguides. Generally, 
good agreement with Eq. (!38l) was found for states with eigenfrequencies / <^ fc- 



B. Coupling to the continuum 

Figure [13] shows the intensity distributions of the third bound state for waveguides with 
bending angles 14.55°, 15.95° and 16.40°. As the bending angle a increases from panel 
(a)-(c), the resonance frequency approaches the cut-off frequency fc and the intensity dis- 
tribution leaks more and more into the leads. This feature is in accordance with Eq. ( I38l) 
since the decay length 5 increases as / approaches fc- For a = 16.40° it was determined 
as 6expt = 50.6 mm from the measured intensity distribution, which is about one quarter of 
the lead length Iq = 211.4 mm. Therefore, it can be expected that the finite lead length 
influences the properties of the state. Indeed, the measured resonance frequency, 9.509 GHz, 
is 8 MHz larger than that of the corresponding waveguide with extension [see Fig. El^c)] the 
length of which is / = 411.4 mm ^ 2lo. Such a shift by several MHz to larger frequencies 
when reducing the lead length was observed for most resonance states close to the cut-off 



frequency and was also predicted in Ref. [27|]. The quality factors Q of the resonances, on the 
other hand, did not change significantly, thus indicating that the leakage due to tunneling 
is small. 

These results demonstrate that for the determination of the eigenfrequencies of the bound 
states the choice of sufficiently long waveguides is crucial. For the experiments presented in 
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Sec. mil the lead length / of the waveguides with extension was chosen significantly larger 
than the decay lengths 6 of the states with frequencies below fc, in order to minimize the 
effect of their finite lengths. 

Next, the transmission of microwave power from port 4 at one lead end through the 



waveguides to port 5 at the other one in Fig. 14(a) was measured in order to investigate 
the coupling of the bound states to the continuum. For this, waveguide-to-coaxial adapters 
(model 4609 by Narda Safety Test Solutions) were attached to the fianges at the ends of the 
waveguides (see Fig. [5]). The adapters close the system but provide a strong coupling to the 



attached coaxial cables. The geometry of this setup is shown in Fig. 14(a) If the coupling 
to the exterior is non- negligible, then the bound states are observable as resonances below 



fc in the measured spectra like in Refs. 
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23 



60|. 



The transmission spectrum obtained by coupling in and out microwave power via the two 



adapters is shown as black thin line in Fig. 14(b) for the waveguide with a = 16.40°, that 
measured between the antennas 2 and 3 located in the leads close to the interior region as 
gray thick line. The latter exhibits three sharp resonances below fc at frequencies which 
were shown to correspond to bound states in Sec. IIIII (see Fig. [H]). The former, on the other 
hand, only shows one resonance below fc, at 9.506 GHz. The intensity distribution (inset 
in the bottom panel) measured at the corresponding frequency confirms that this resonance 
indeed can be assigned to the third bound state which apparently couples to the adapters. 
This can be attributed to its large decay length in the leads. Thus the coupling of the third 
bound state to the continuum indeed implies a resonance in the transmission spectrum due 
to resonant tunneling, whereas the first and second bound state do not. This is explainable 
with their considerably smaller decay lengths. 

Note that the resonance frequency of the third bound state differs slightly for the mea- 
surements with the antennas and the adapters, respectively. This is attributed to the fact 
that the coupling of the adapters to the waveguide is much stronger than that of the small 
wire antennas, as can be deduced from the average value of jS'^ap above fc, which is larger 
in the former case. 
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C. Summary 



Figure 14(b) shows that the bound states can indeed modify the transmission spectrum 



of the bent waveguides (see also Ref. |59|) due to resonant tunnehng as observed in quantum 



wires [22|]. However, this was seen in just two cases for all nine waveguides, and only if no 
extensions for the leads were used. In both cases, the associated resonance frequency was 
less than 5 MHz apart from the cut-off frequency. Apparently, for resonant tunneling to 
be observed, the decay length 6 of the corresponding state must be about the same as the 
lead length. However, only states close to the cut-off frequency in waveguides with bending 
angles close to a critical angle satisfy this condition which emphasizes the importance of the 
knowledge of the critical angles. 

Furthermore our experiments show that especially for these cases the finite lead length 
and the details of the coupling can become important, which requires further theoretical 
investigations. First results regarding the influence of the lead length were published in 
Ref. y. 

VI. CONCLUSIONS 

We have studied experimentally and numerically the occurrence of bound states in sharply 
bent waveguides below the threshold for propagating modes. In some of the experimental 
realizations the bending angle was chosen close to a critical angle. The measured eigen- 
frequencies of the bound states and the values of the critical angles were predicted to a 
high accuracy with two independent methods: The first method consists of the numerical 
solution of the Schrodinger equation for a free particle with Dirichlet boundary conditions 
at the walls of the waveguide based on a scattering approach. 

The second one encompasses the derivation of a Hamiltonian containing a binding po- 
tential by introducing a suitable conformal mapping. The associated Schrodinger equation 



was solved semiclassically using corner-corrected WKB approach [3l|. These results are 
of particular relevance for the transmittance of electromagnetic waveguides and quantum 
wires. The theoretical treatment is, to our knowledge, the first one to work out the problem 
by analytical methods, even though it relies on an approximation. The demonstration of its 
applicability to sharply bent waveguides renders possible the analytical treatment of two- 
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dimensional problems with non-trivial boundaries through the use of a conformal mapping 
and an extended WKB approximation. 

In further experiments, the intensity distributions of the bound states were measured. 
They exhibit interesting features close to the cut-off frequency and for bending angles close 
to a critical one. Especially the influence of the finite lead lengths on their boundedness, 
that is, the size of their coupling to the exterior, was investigated. These studies revealed 
that predominantly bound states with resonance frequencies close to the cut-off frequency 
are affected if the length of the leads is of the same size as their decay length, the relevant 
mechanism being resonant tunneling. 
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Appendix A: Numerical computation of bound states 

The Schrodinger equation corresponds to that of a free particle with Dirichlet boundary 
conditions along the boundary dQ, 

(A,,y + A;2)0(x,i/) = 0, 0(x,i/)|an = O, (Al) 

where k is the wave number related to the energy E in Eq. ([T]) as E = k"^. The waveguide 
consists of two leads 1 and 2 with parallel walls and an interior region, as depicted in Fig. [151 
The coordinate system (x, y) is chosen such that its origin is located at the outer corner, 
and the x-axis coincides with the symmetry axis of the bent waveguide. In a scattering 
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experiment waves are sent into the interior region via one lead and exit either via the same 



)in g th is process is determined with the 
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39 



40 



44|- With that method, 



or via the other one. The scattering matrix descri 
help of the scattering formalism applied in Refs. 
the wave functions y) solving Eq. flAip are obtained by writing down general solutions 
valid in the leads and in the interior region, respectively, and joining them smoothly along 
the common chord. 

The solutions (j){x, y) of Eq. flAip are either symmetric or antisymmetric with respect to 
the X-axis, that is, either (p{x,y) = (j){x,—y) or (p{x,y) = — — y). Hence, they vanish 
in the antisymmetric case along the x-axis, (j){x, 0) = 0. Thus y) obeys the Dirichlet 
boundary condition not only along the walls of the waveguide but also along the chord y = 0. 
This fact defines a further wall which divides the bent waveguide into two straight ones with, 
respectively, only one open end. Consequently, waves entering one of its leads are eventually 
reflected specularly at the end wall and then exit again, i.e., all antisymmetric states are 
unbound. Since our aim is the computation of the wave numbers and the wave functions of 
the bound states, which have been shown to exist below the cut-off frequency of the first 
propagating mode, we consider only the symmetric solutions of Eq. ( lAll) . 

We introduce in each lead a coordinate system with origin at the inner corner of the 
bent waveguide and one axis along the inner wall of the respective lead and the other one 
perpendicular to it (see Fig. [15]). Let L = 1 denote the lead with y > in Fig. [15] and L = 2 
that with y < 0, W the width of the leads and d = Wcsc{a/2) the distance between the 
outer and the inner corner of the bent waveguide. 



With respect to the coordinate system {x,y) the coordinates {xi^yL) are given as 



XL 



cos (fL SmifL 

^ — sin (fL cos ifL 




(A2) 



where ^Pl=i = «/2 and ^Pl=2 = 37r/2 — a/2. 
We introduce the notation 

k - 



/ mrr 



(A3) 



where m = 1,2,... denotes the transverse mode numbers in the leads. The wave number 
km is purely imaginary for A; < ^ and thus does not correspond to a propagating mode. 
For a given k the number Mo of the latter equals the integer part of Wk/n, i.e., 

'Wk' 



TT 



(A4) 
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Below the cut-off wave number kc of the first propagating mode, i.e., for k < kc = no 
propagating modes exist in the leads. 

The solutions of Eq. (1A1I) in the leads are given by 



Vlk^l \/w/2 



(A5) 



for lead 1, i.e., y > 0, and 



(A6) 



for lead 2, i.e., ?/ < 0. Here, Li, L2 G {1,2} refer to the leads, and 5l2M — -^1 7^ -^2, 
^L'zM ~ ^ -^1 ~ -^2 and 6'(x) = 1 for x > and zero otherwise. 

In both equations the first term corresponds to an incoming wave with mode number 
n propagating from the asymptotic region, i.e., the ends of the leads, into the interior 
region, while in the second term the sum runs over all Mo outgoing propagating and 
exponentially decaying modes with a purely imaginary km- Since these evanescent modes 
couple only weakly to the interior region for large mode numbers m, their total number 
can be limited to some finite value M^. The reduced T-matrix with n,m < Mg yields the 
scattering matrix S^^'^^. The matrix element 5*^^^^ describes the scattering from mode n 
in lead Li via the interior region into mode m of lead L2. 

In the interior region, the solutions of Eq. ( lAll) . which are symmetric with respect to the 
X-axis and vanish along its outer boundary, that is, along its common boundary with the 
bent waveguide, can be written explicitly in polar coordinates 

X = rcosO 

(A7) 

y = rsin6' 

as 

Mt 

"^n'ir, e) = J2 B^.n COS m Ji,{kr), (A8) 

where we introduced the notation /t = ^(2/i — 1). Due to the property of the Bessel functions 
Jfi{x) that they become vanishingly small for fi ^ x, the infinite expansion can be truncated 
at a certain value n = Mt, which depends on that of kr. 
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The matrix elements T^^^^ entering the ansatz for \ ^^"^ ^i^L inside the lead region 



given by Eqs. (]A5|) and (]A6P and the expansion coefficients -B^„ entering \E'™*(r, ^) for the 
interior region in Eq. (lASP are chosen such that the wave functions \E'^*(r, 6) and \l'^^^(a;i, 
respectively, \E'^ j;^(x2, 1/2) and their normal derivatives match at their common boundaries 
represented in Fig. [15] by the dotted lines and defined by (xi = 0, < ?/i < W) and 
(0 < X2 < W, 1/2 = 0), giving rise to the conditions 



^i-*(xi = 0,1/1) = ^1^1^ (xi = 0,yi) 
^int(a;2,t/2 = 0) = <'^(a;2, 1/2 = 0) 



(A9) 
(AlO) 



and 



(1) 



9^2 



-(3^2,2/2) 



xi=0 



S/2=0 



(2) 



dy2 



xi=0 



1/2=0 



(All) 
(A12) 



Here, \E'™*(xl, ?/l) is expressed in terms of the coordinates of lead L. To compute the partial 
derivatives of ^'™*(x,|/) defined in Eq. flASP with respect to these coordinates we use the 
relations 



a 

dxi 



cos(a/2-^)|: + 



d I sin(a/2-6>) d 



dr 

a 



sin(a/2+e) a 
r (96»' 



(A13) 



^ = cos(a/2 + ^^)f - 
with the polar coordinates r, 6 defined in Eq. (IA7p . 

Multiplying both sides of Eq. f lATT]) with sin (^l/i) and of Eq. flAT2D with 

sin (^^X2^ I \JW/2 and integrating them along the common chord leaves exactly one matrix 
element on the right hand side of the respective equation [see Eqs. (1A5I) and (1A6I) ]. We 
introduce the (Mo + M^) x Mt dimensional matrix A with elements 



A 



v/|fcn| ra/2 1^ Wcot(a/2) 
ifcn Jo '^'^cos2 (a/2-61 ; 



(a/2-e) 



W/2 



. w 



X 



cos(/i^) cos(a/2 - ^)A; 



x=kp 



(A14) 



sin(/i^) sin(a/2 - ey^^J~^{kp) 



where 



with 



^ = -W cot(a/2) tan(a/2 -e) + W , 



iycot(a/2) 
cos(a/2 - 9) 



(A15) 
(A16) 
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forn = 1, Mo + Mc and /i = 1, Mt. 

Then we obtain for Li, L2 G {1, 2} the matrix equations 

- 15l,m + =A-B. 



(A17) 



They relate the matrix elements T^^^^ of T^^'^^ to the matrix elements B^n of B. The 
values of the latter are fixed by a further set of matrix equations. To derive it, we introduce 
the arc length s, which varies along the perimeter of the interior regiori, and expand the 
wave function \E'™*(x, with x, y G dVt in a Fourier series in s (cf. Ref. |40|). Here, s = yi 
{s = —X2) along the common boundary with the lead L = 1 (L = 2), and s = r along 
that with the wall of the bent waveguide at 6 = ±a/2 in terms of the polar coordinates 
defined in Eq. flA7p . Using the fact that \E'™*(r, 6) vanishes along the outer walls, the Fourier 
coefficients associated with s read 



D 



"/2 W^cot(a/2^ /27rA/ 
^ ^os^(a/2-g) ^^^^^^ I 



r 



(AI8) 



with A, /i = 1, . . . , Mf. The variables ^ and p are defined in Eq. f lA16P and F = 214^ [1 + 
cot (a/2)] denoting the length of the perimeter of the interior region. Defining the Mt x 
(Mo + Mc)-dimensional matrix C with elements 



w 1 



/in 



l—cos[-iTn+2iTii- 



1— cos(7r?i— 27r/i-^ 



with n 



7rn— 27r/i p 

Mf, Eqs. dSni) and " 



Mo + Mo, /i = 1, . 



(A19) 



yield the matrix equations 



2D-B, 



(A20) 



for L e {1,2}. 

Combining Eqs. (IA17P and flA20p . the matrix elements of B are obtained as 



B= D-CA 



C. 



(A21) 



For k > kc, that is, if the leads allow for propagating modes, we obtain with Eq. (1A16P 
setting Mo = the final result for the scattering matrices 



^1,1 
^1,2 



^2,2 
^2,1 



1 + AiD-CA 

-1 



-1 



C 



i (z) - cA) c. 



(A22) 
(A23) 
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The full scattering matrix 

^^^•1 S 



( cl,2 



(A24) 



V 



^2,1 ^2,2 

is symmetric according to Eq. f lA23p . Using that C and D are real and that A is purely 
imaginary, it can be shown that 5* is unitary. 

Below the cut-off frequency, that is, for k < kc the first term on the right hand side of 
Eqs. flA5l) and flA6l) and on the left hand side of Eqs. flA17|) and flA20l) vanishes, leading us 
to the equation 

(l) - CAj B = (A25) 

for B, or equivalently 

det (z) - CA) = 0. (A26) 

The matrix elements of A, C and D depend on k. Those values of k which solve Eq. flA26P 
correspond to the wave numbers of the states, which are bound in the interior region. The 
associated wave functions are obtained by inserting these values into Eqs. flASp and flA6p . 
omitting the first term in these equations, and into Eq. (lASp . In these calculations we 
typically chose for the number of closed channels Mc and the number of terms Mt in the 
sum Eq. flASP values between 1 and 5. 
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(a) (b) 
FIG. 1: Analogy between the decay of a compound nucleus (a) and a bound state in a bent 



waveguide |(b)[ A compound nucleus is formed by bombarding a nucleus A with a particle a 
and eventually decays into a residual nucleus and particles, B + h, C + c, etc.. This coupling to 
the exterior region is modeled by attaching different formation and decay channels to a confined, 
interior region. Panel |(b)| shows a sketch of a quantum wire or electromagnetic waveguide with a 
sharp bend of angle a, the straight parts of which have a length / and a width W. The border 
between them and the corner region is indicated by dashed lines. Since bound states form in the 
latter, we also denote it as interaction region. The coordinate along the leads is denoted by x^. 
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FIG. 2: Numerical results for the rescaled eigenfrequencies ///c with / the eigenfrequencies of the 



bound states and fc the cut-off frequency versus the bending angle a of the waveguide. Figure 2{a] 



shows the ratios for 0.2° < a < 4°, Fig. [2(b)] those for 4° < a < 60°. Each curve corresponds to the 
eigenfrequency of one bound state in units of fc, which is a continuous function of the bending angle 
a, starts at a = 0° and increases until it reaches the value ///c = 1 at a certain bending angle. 
There the bound state turns into a propagating mode. The corresponding angles are called critical 
angles acrit- The arrows in panel (b) indicate the critical angles considered in the experiments with 
microwave waveguides. The number of bound states increases rapidly as a approaches 0°. Indeed, 
for a w 0.2° the curves are barely distinguishable. The eigenfrequency of the first bound state, 
which corresponds to the lowest curve in both panels, converges to f / fc = 1/2. This is the cut-off 
frequency of a waveguide of double the width of the leads. The dots mark the eigenfrequencies of 
the four bound states at a = 9.65°, whose intensity distributions are shown in Fig. [3l 
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FIG. 3: (Color online) Square of the numerically computed wave functions associated with the four 



bound states at a = 9.65°. The corresponding eigenfrequencies are marked by dots in Fig. 2(b) 



Panel (a) shows the intensity distribution for ///c = 0.631, panel (b) that for ///c = 0.748, panel 
(c) that for ///c = 0.851 and panel (d) that for ///c = 0.946. The blue (darkest) color corresponds 
to the lowest, the yellow (brightest) to the highest intensity. 
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FIG. 4: Number of bound states N{a) in its dependence on the bending angle a. The circles 
show the numerical results, the solid line is obtained from the analytical expression Eq. The 
agreement between the numerical results and the analytical expression is very good. The inset 
shows a magnification for Ocrit < 1°- 
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FIG. 5: (Color online) Photo of a bent waveguide with angle a = 26.50°. The three small copper 
blocks soldered on top of the waveguide are used to attach the antenna ports (not shown) labeled by 
numbers 1, 2 and 3, where microwave power is coupled into and out of the waveguide. Additional 
straight waveguide parts can be attached to the two flanges. 
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FIG. 6: Bound states of a bent waveguide, here with bending angle a = 16.40°, identified by the 
measured resonance spectrum shown in logarithmic scale. The gray bar indicates the frequency 
range where the cut-off frequency fc is expected. The broad resonances observed above fc are 
due to partial reflections of the microwaves at the open ends of the waveguide. The intensity 
distributions for the resonances labeled (a)-(d) are shown in the corresponding panels of Fig. [71 
The arrow indicates the position of the resonance (d), which is above the cut-off frequency fc and 
thus broad. The inset is a magnification around fc including resonance (c). 
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FIG. 7: (Color online) Intensity distributions of three bound states and one propagating mode for 
the waveguide with a = 16.40° and length / = 411.4 mm measured with the perturbation body 
method. The electric field intensity / oc l-E^P is plotted in false colors. The blue (darkest) color 
corresponds to the lowest, the yellow (brightest) to the highest intensity. Panels (a)-(d) correspond 
to the resonances indicated in Fig. [HI In panels (a)-(c), i.e., for / < /c, the intensity distribution 
is localized around the inner corner, whereas for / > /c it periodically extends over the whole 
waveguide as shown for one example in panel (d). In panels (a) to (c) yellow (bright) stripes are 
observed, and the intensity does not seem to vanish at the walls of the waveguide. This is an 
experimental artifact due to a lag of the perturbation body behind the guiding magnet during the 
measurements. Furthermore, the intensity distributions do not show the expected symmetry. This 
again is an experimental artifact resulting from the choice of the rectangular grid of positions of 
the perturbation body, one side being parallel to, respectively, the lower lead of the waveguides in 
the figure. 
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FIG. 8: Measured (crosses) eigenfrequencies of the bound states in units of the cut-off frequency 
in their dependence on the bending angle a. The sohd hues correspond to the eigenfrequencies 



computed in Sec. |TI] with a conventional numerical method (same data as in Fig. 2(b), but shown 
only in the range 8° < a < 30°). When an eigenfrequency crosses the dashed line, it reaches the 
value of the cut-off frequency. The occasional, slight deviations between the experimental and the 
numerical results are attributed to imperfections in the construction of the waveguides resulting 
from the soldering process. Also small differences between the critical angles computed in Sec. HI] 
and the experimental ones are observed. 
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FIG. 9: Grid of the conformal coordinates {x{u,v),y{u,v)) obtained from the map z = -F(^) given 
in Eq. ([8]) for a = 180°, 126° and 54°. For a fixed v = const, the coordinate lines are perpendicular 
on the boundaries (x(n = 1, v),y{u = 1, v)) and {x{u = 0, v), y{u = 0, v)), while for fixed u = const, 
the coordinate lines are parallel to the latter for large values of v. An abrupt change of the contour 
density occurs near the corners which is related to the Jacobian of the transformation z = -F(^) as 
expressed by Eq. (I15p . This feature stands outmost clearly in the example shown on the right with 
an accumulation of lines at the inner corner and a blank space to the right of the outer corner. 
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FIG. 10: [(a)] Effective potential K=i(v) = 1 - T4=i(w) and [(b)] the derivative ^l^i^iM 



dv 



dv 



for Q = 1° (dotted line), a = 25° (dashed-dotted line), a = 49° (dashed-dotted-dotted line), 
a = 73° (dashed-dashed-dotted line), a = 97° (dashed line) and a = 121° (full line). 
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FIG. 11: (a) Comparison between the rescaled eigenfrequencies f / fc of the bound states computed 
with the conventional numerical method outlined in Sec. |TI] (fuU hues) with those obtained from 
the corner-corrected and ground-state shifted WKB approximation (dots). |(b)| Relative deviations 
of the critical angles evaluated numerically and with the WKB approximation versus a^^^ . Both 
panels demonstrate that the agreement between the numerical and the WKB results is already 
good for the third bound state. 
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FIG. 12: (Color online) Decay of a bound state into the leads represented by the projection of the 
measured field distribution (gray line) onto the x^-axis parallel to the waveguides [see Fig. l(b)| . 
Here we use the second bound state at 9.364 GHz for a = 27.65° with lead length /q and show this 
projection in linear scale (top panel) and logarithmic scale (bottom panel). In the regime from 
20mm < xj, < 100mm the intensity distribution decays exponentially, indicated by the dashed 
line, with the decay length (5expt = (16.2 it 0.4) mm. The inset displays the complete intensity 
distribution but parts of the waveguides are not shown. 
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FIG. 13: (Color online) Influence of the short leads on the bound states as expressed by the 
leakage of the measured intensity distributions of the third bound state into the leads for angles 
(a) a = 14.55°, (b) a = 15.95° and (c) a = 16.40°. The electric field intensity / oc \E^\'^ is 
plotted in false colors. The blue (darkest) color corresponds to the lowest, the yellow (brightest) 
to the highest intensity. In contrast to Fig. [71 the waveguides now have short leads with length 
Iq summarized in Tab. HI With increasing bending angle the intensity penetrates deeper into the 
leads. 
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FIG. 14: (Color online) Resonant tunneling in a bent waveguide observed in measured microwave 



transmission spectra, (a) Sketch of the setup. The crosses marked 2 and 3 indicate the positions 
of the wire antennas also used in the experiments presented in Sec. IIIII and the crosses marked 4 
and 5 denote the positions of the waveguide-to-coaxial adapters. |(b)| Frequency spectrum of the 
waveguide with a = 16.40° and adapters attached. The spectrum was measured both with the 
adapters (black thin line) and with antennas 2 and 3 located near the tip of the waveguide (thick 
gray line). The light gray vertical bars and its width indicate the cut-off frequency fc and its error 
bar, respectively. The two resonances at 6.49 GHz and 8.15 GHz observed in |532p are not seen in 
1 5*54 P because the corresponding intensity distributions decay rapidly in the leads and thus cannot 
couple to the adapters. Only the third resonance below fc at 9.506 GHz is observed in 1 5*541^, since 
its intensity distribution extends far into the waveguide (see inset in the bottom panel). 



47 




FIG. 15: Definition of the two coordinate systems used in the analysis of a sharply bent waveguide. 
The origin of the coordinate system (x, y) coincides with the outer corner and the x-axis with the 
symmetry line of the waveguide. Moreover, we identify two asymptotic regions defined by the leads 
and one interior region around the inner and the outer corner of the waveguide which suggests to 
introduce two coordinate systems (xi,yi) and (x2,y2) defined in the leads L = 1 and L = 2, 
respectively. 
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Tables 
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TABLE I: List of the bending angles a of the waveguides investigated experimentally. The lengths 
/ and /() (^f t he; leads are with and without an additional waveguide extension, respectively. 



a (deg) 


/q (mm) 


/ (mm) 


28.40 


122.2 


422.2 


27.65 


356.7 


456.7 


26.50 


131.6 


431.6 


16.40 


211.4 


411.4 


15.95 


224.8 


424.8 


14.55 


240.0 


440.0 


11.83 


300.1 


400.1 


11.00 


328.2 


428.2 


9.65 


361.2 


461.2 
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